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Abstract 



We have searched for direct production of scalar top quarks at the Collider 
Detector at Fermilab in 88 pb _1 of pp collisions at yfs =1.8 TeV. We assume 
the scalar top quark decays into either a bottom quark and a chargino or a 
bottom quark, a lepton, and a scalar neutrino. The event signature for both 
decay scenarios is a lepton, missing transverse energy, and at least two fr-quark 
jets. For a chargino mass of 90 GeV/c 2 and scalar neutrino masses of at least 
40 GeV/c 2 , we find no evidence for scalar top production and present upper 
limits on the production cross section in both decay scenarios. 
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The Minimal Supersymmetric extension to the Standard Model (MSSM) [Ij assigns a 
scalar supersymmetric partner for every Standard Model fermion and a fermionic superpart- 
ner for every Standard Model boson. The weak eigenstates of each scalar superpartner mix, 
forming mass eigenstates J|. The splitting of the mass eigenvalues is proportional to the 
mass of the Standard Model partner. Therefore, the superpartners of the top quark weak 
eigenstates, It and Ir, may have the largest mass splitting of all the scalar quarks (squarks). 
The running of the squark mass parameters is proportional to the Yukawa coupling of the 
Standard Model partners, such that the diagonal elements of the tx,tn mass matrix should 
be smaller than those of the other squarks 0. Thus, the lighter scalar top mass eigenstate, 
ti, is the best candidate for the lightest squark and is potentially lighter than the top quark. 
We report the results of a search for direct production of tit± in 88±4 pb _1 of data collected 
during the 1994-95 Tevatron run using the Collider Detector at Fermilab (CDF). 

The CDF detector has been described elsewhere ||. In this analysis, we used electrons 
identified in the central electromagnetic calorimeter which covers the pseudorapidity region 
1 77 1 <1.1. We used muons identified by tracks in drift chambers in two detector subcompo- 
nents outside the calorimeters. The first muon subsystem is located behind five absorption 
lengths of material and covers the region \rj\ <0.6. The second is located behind an additional 
three absorption lengths of material and has the same rj coverage as the first. 

Scalar tops are strongly produced in the Tevatron via qq annihilation and gluon-gluon 
fusion. We searched for t\ti production within the framework of the MSSM for the case where 
mi x < m t . We assumed R-parity [0] is conserved and restricted ourselves to two separate 
t\ decay modes Q. In the first, the decayf] t\ — > bxl , where xt is the lightest chargino, 
proceeds with a branching ratio of 100%. We required one of the charginos, which decay 
via a virtual W, to decay as xf ~~ > e+h> Xi or / i+z/ Xi> where Xi is the lightest neutralino, 
with an assumed branching ratio of 11% for each lepton type. For models where t% — > bxf is 
not kinematically allowed, we considered a second decay scenario in which t\ —>■ bl + i), where 
v is a scalar neutrino and each I — e, fi, r has a branching ratio of 33.3%. In these two 
scenarios, either the x° or the v is the lightest supersymmetric particle (LSP) and does not 
decay. A third possible decay scenario in which the t\ — > cx\ branching ratio is 100% is the 
subject of separate CDF searches 0. 

In both decay scenarios considered here, the t\ signature is at least one isolated lepton, 
missing transverse energy (j^r ) from the neutral LSP's and at least two jets from the b 
quarks. This signature is very similar to that of the top quark, with kinematic differences 
due to the smaller t\ mass in our search region, the presence of two massive neutralinos in 
the final state, and the absence of a real W in the final state. We therefore expect events 
with lower lepton px, lower jet and multiplicity, and without a peak in the lepton- 
$t transverse mass. To remain efficient for the smaller t\ mass, we used data collected with 
the low-pr electron and muon triggers described in Ref. ||. These trigger thresholds were 
Et > 8 GeV for electrons and px > 8 GeV/c for muons. 

The data for this analysis were obtained by requiring (i) an electron with Et > 10 GeV 
or muon with pt > 10 GeV/c originating from the primary vertex and passing lepton 
identification cuts, (ii) St > 25 GeV, and (iii) at least two jets with cone sizes of 



Unless otherwise noted, decay channels imply their charge conjugates. 
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R = ^J{Ar]) 2 + (A</>) 2 = 0.7, one with E T > 12 GeV and the second with E T > 8 GeV. 
The lepton identification cuts were identical to those used in previous CDF analyses |7||§. 
For electron identification, the electron was required to have lateral and longitudinal shower 
profiles consistent with those of an electron, have less than 5% of its energy deposited in the 
hadronic calorimeter, and be well-matched to a track from the CTC. A muon was required 
to have tracks in the inner and outer central muon chambers which were well-matched to a 
track from the CTC. We further required the leptons to pass an isolation cut in which the 
calorimeter Et in a cone of R = 0.4 around the lepton was less than 2 GeV (excluding the 
lepton tower). No explicit tau identification was conducted, but leptons from tau decays 
were accepted. 

We used the SVX 7 detector to identify secondary vertices from b quark decays and selected 
events with at least one secondary vertex. The tagging algorithm is described in Ref. 



with improvements given in Ref. || and efficiency measured in Ref. |10|. We reduced the 
Drell-Yan background in our sample by removing events with two isolated, opposite-sign 
leptons. This background was further reduced by removing events with an isolated lepton 
that reconstructed an invariant mass > 50 GeV/c 2 with any additional, isolated CTC track. 
Finally, we reduced the background from bb events and events with hadrons misidentified as 
leptons (fake leptons) by requiring that the A(j) between the $t and the nearer of the two 
highest- Er jets be > 0.5 rad. This reduces fake $t due to jet energy mismeasurement. The 
number of events remaining in our sample after all cuts is 81. 

Signal and background selection cut efficiencies were estimated using a variety of Monte 
Carlo generators followed by a CDF detector simulation. Signal event samples were created 
using ISAJET version 7.20 ||TT|. The supersymmetric particle masses used in signal simulation 
were: m-± = 90 GeV/c 2 , m^o = 40 GeV/c 2 , and rric, > 40 GeV/c 2 , which are consistent 
with current lower limits ||12|| . The signal selection efficiency increases with mi but decreases 
with m,ji and m-±, reaching a plateau as event energies advance from cut thresholds. Some 

specific efficiencies are 5.4% for t\ — > bl + v (m^ = 130 GeV/c 2 , m„ = 40 GeV/c 2 ) and 0.7% 
for ti -> bxt {m h = 120 GeV/c 2 , m-± = 90 GeV/c 2 , and m^o = 40 GeV/c 2 ). These 
selection efficiencies include branching ratios of forced decays. 

The significant sources of uncertainty for signal selection efficiency are (i) the 6-jet tagging 
efficiency, (ii) the trigger efficiencies, (iii) the luminosity, and (iv) initial- and final-state 
radiation. The effects of some of these sources vary with , but none contribute more 
than 10% to the overall uncertainty, which is less than 16% for all m^ considered. 

Standard Model backgrounds come from any process that can produce two or more 
jets, either real or fake leptons, and real or fake $t. This includes heavy flavor quark 
production, vector boson production with two or more accompanying jets, and inclusive jet 
production with real or fake leptons. The number of events from the first two processes that 
we expected in our data sample were predicted using measured or calculated cross sections 
and selection efficiencies determined from Monte Carlo. Top-pair and single-top production 



were simulated using HERWIG version 5.6 |13[]. For m t = 175 GeV/c a t i is 5.1 ± 1.6 pb 



and a t i for VF-gluon fusion from a next-to-leading-order (NLO) calculation is 1.70 ± 0.15 pb 



14| . Vector boson samples were generated using VECBOS version 3.03 [1J| and normalized 
according to CDF measurement [|T^]. Drell-Yan, bb, and cc samples were generated with 
ISAJET version 7.06 and normalized to independent CDF data samples. 

To determine the number of events with fake leptons in our sample, we used a data 
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sample passing all our selection cuts with the exceptions of a non-overlapping Et requirement 
(15 < Et < 20 GeV) and no requirement on A(/)(Et, nearer jet). The number of fake lepton 
events was normalized to this data sample, which contained negligible signal, after other 
backgrounds were subtracted. The number of fake lepton events was then extrapolated to 
the signal region using cut efficiencies determined from an independent fake-lepton event 
sample. 

The complete list of backgrounds and the number of expected events remaining after all 
cuts is given in Table |. The significant backgrounds are it, bb, W ± (—> l ± i')+ >2 jets, and 
fake lepton events. The number of data events agrees well with the expected background. 

To determine the number of potential signal events in this final data sample, we per- 
formed extended, unbinned likelihood fits for each t\ mass considered for both decay scenar- 
ios. The likelihood fits compared the shapes of distributions of the signal and background 
and included Gaussian terms tying the fit background levels to their predicted levels. The 
fit parameters were the numbers of signal events, ft events, bb plus fake lepton events, 
and vector boson events (represented in the fit by the W ± + >2 jets distributions). We 
used the Kolmogorov statistic applied to the simulated distributions of signal and combined 
backgrounds to determine the most sensitive kinematic distributions to use in the fit. The 
kinematic distributions evaluated include lepton p T , H T (the scalar sum of lepton E T , Et , 
and jet E T for all jets with E T > 8 GeV), jet multiplicity, and A0(jetl,jet2), where jets are 
ordered in E T . 

For the t\ — > bxt decay, sensitivity to signal was greatest for a two-dimensional fit to the 
combined probability distributions for Ht and A0(jetl,jet2). Fit results at all masses were 
consistent with zero signal events. The fit result for t\ — > bxt with = 115 GeV/c 2 , m-± 
= 90 GeV/c 2 , and m^o = 40 GeV/c 2 is shown in Fig. 0. The 95% C.L. limits on a- t ~ h for 
this decay are shown in Fig. |^ as a function of [|T7|]. The NLO prediction for cr^- using 
the renormalization scale fi = and parton distribution function CTEQ3M is shown in 
Fig. ^| for comparison [Pf|. 

For the ti — > bl + v decay scenario, sensitivity to signal was greatest for a fit to the Ht 
distribution. Again, all fit results were consistent with zero signal events. The 95% C.L. 
limits on a^j for the t\ — > bl + u decay are shown in Fig. [3] for = 40 and 50 GeV/c 2 . 
We consider the regions of supersymmetric parameter space for which the 95% C.L. limit 
on is less than the NLO prediction (/i = rajj to be excluded. The resulting excluded 
region in the plane of mi versus m„ for is shown in Fig. [|. 

To conclude, we have searched for direct titi production in 88±4 pb _1 of data collected 
using the CDF detector during the 1994-95 Tevatron run. We found no evidence for titi 
production for either t\ — > bxt or ti ~^ bl + v and present upper limits on = as a function 
of m h . 
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Council of Canada; the National Science Council of the Republic of China; the A. P. Sloan 
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TABLES 



TABLE I. Number of data events and expected background events after all selection cuts. 
The dominant sources of uncertainty on the numbers of expected events are integrated luminosity, 
cross section, trigger efficiency, and b-jet tagging. (Fake leptons are hadron tracks which have been 
misidentified as leptons.) 



process 


number of events 
expected after all cuts 


W ± (-^ e ± v or n ± v)+ >2 jets 


44.5 ± 7.3 


tt 


17.8 ±4.5 


bb 


5.8 ±0.8 


W ± (-^ t ± u)+ >2 jets 


2.6 ±0.4 


tb (from W — g fusion) 


1.6 ±0.2 


Z(-> e + e~ or n + n~)+ > 2 jets 


1.4 ±0.2 


Z{-+ t + t-)+ > 1 jet 


0.4 ±0.1 


7 -► i+r 


0.4 ±0.1 


cc 


0.06 ± 0.02 


fake lepton events 


12.7 ± 1.6 


background total 


87.3 ±8.8 


data 


81 
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FIG. 1. Results of the two-dimensional fit to and A</>(jetl,jet2) when the t\ — > b\i 
branching ratio is 100%, m fl = 115 GeV/c 2 , m~± = 90 GeV/c 2 , and m^o = 40 GeV/c 2 . The 
quantities Ht and A^>(jetl,jet2) are defined in the text. The points represent the data. Cumulative 
contributions from bb and fake lepton events, tt, and — > Z 1 * 1 ^ + jets are represented by dotted, 
dashed and solid lines respectively. There is no significant contribution from signal. To illustrate 
the shape difference, a signal distribution with arbitrary normalization has been overlaid with a 
dot-dash line. 
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FIG. 2. The points represent the CDF 95% C.L. cross section limit as a function of t\ mass 
when the t\ — > b\t branching ratio is 100%, m~± = 90 GeV/c 2 , and m^o = 40 GeV/c 2 . The line 
without markers represents the NLO prediction for o~ ~ using the renormalization scale \i = . 
The dashed lines represent the NLO cross section for \i = and jjl = 2m^. 
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FIG. 3. CDF 95% C.L. cross section limit as a function of t\ mass when the t\ — > 
branching ratio is 100% and mj; = 40 GeV/c 2 (squares) or 50 GeV/c 2 (triangles). The line 
without markers represents the NLO prediction for o~ ~ using the renormalization scale /U = . 
The dashed lines represent the NLO cross section for \i = and jjl = 2m^. 
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FIG. 4. 95% C.L. excluded region in the plane of versus rriy when the t\ — > be + u, 
t\ — > bn + u, and t% — > bT + v branching ratios are 33.3%. We define the exclusion region as that 
region of supersymmetric parameter space for which the 95% C.L. limit on | is less than the 
NLO prediction (/x = m^). The LEP1 rriy limit and ALEPH excluded region in the m t - versus 
m u plane are also shown p9[ . The ALEPH excluded region corresponds to the case in which the 
ti decouples from the Z°. 



13 



